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Abstract. We investigate the evolution of the 0.5-2 keV 
soft X-ray luminosity function (SXLF) of active galac- 
tic nuclei (AGN) using results from ROSAT surveys of 
various depth. The large dynamic range of the combined 
sample, from shallow large-area ROSAT All-Sky Survey 
(RASS)-based samples to the deepest pointed observation 
on the Lockman Hole, enabled us to trace the behavior 
of the SXLF. The combined sample includes about 690 
AGNs. As previously found, the SXLF evolves rapidly as 
a function of redshift up to z ^ 1.5 and is consistent with 
remaining constant beyond this redshift. 

We have tried to find a simple analytical description 
of the SXLF in the overall redshift and luminosity range, 
using Maximum-Likelihood fits and Kolgomorov-Smirnov 
tests. We found that a form of the Luminosity-Dependent 
Density Evolution (LDDE), rather than the classical Pure 
Luminosity Evolution (PLE) or the Pure Density Evo- 
lution (PDF) models, gives an excellent fit to the data. 
Extrapolating one form of the LDDE model (LDDEl) ex- 
plains « 60% of the estimated soft extragalactic Cosmic 
X-ray Background (CXRB). We have also found another 
representation (LDDE2), which produces « 90% of the 
CXRB and still gives an excellent fit to the sample AGNs. 
These two versions of the LDDE models can be considered 
two extremes of the possible extrapolations of the SXLF 
below the flux limit of the survey. 

We have also investigated the evolution of the number 
density of luminous QSOs with LogL^ > 44.5 [/i5oergs~^], 
where the evolution can be traced up to the high red- 
shift. We have compared the results with similar quanti- 
ties in optically- and radio-selected luminous QSOs. Un- 
like these QSOs, evolution of the i?OS'^T'-selected QSOs 
do not show evidence for the decrease of the number den- 
sity in z ^ 3. The statistical significance of the difference 
is, however, marginal. 
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1. Introduction 

The AGN/QSO luminosity function and its evolution with 
cosmic time are key observational quantities on under- 
standing the origin of and accretion history onto super- 
massive balckholes, which are now believed to occupy the 
centers of most galaxies. Since X-ray emission is one of 
the prominent characters of the AGN activity. X-ray sur- 
veys are effective means of sampling AGNs for the lumi- 
nosity function and evolution studies. The Rontgen satel- 
lite (ROSAT), with its unprecedented imaging capabil- 
ities, provided us with soft X-ray surveys with various 
depths, ranging from the i?05^T All-Sky Survey (RASS) 
to the i^O^^r Deep S urvey (RDS) on the Lockman Hole 
(Hasinger et al. 1998 ). Various optical identification pro- 
grams of the survey fields have been conducted and the 
combination of these now enabled us to construct the soft 
X-ray luminosity function (SXLF) as a function of red- 
shift. 

The evolution of SXLF has already been seen in 
the Extended Me dium Sensitivity Survey (EM SS) AGNs 
(Maccacaro et al. 1991; Delia Ceca et al. 1992| ) for high- 
luminosity AGNs. Combining results from deep ROSAT 
PSPC surveys and the EMSS has extended the sample 
into the higher-redshift lower-luminosity regime, provid- 
ing much wider baseline to explore the evolution proper- 
ties (e.g. Boyle et al. 1994; Jones et al. 1996; Page et al. 
1996). All of these were characterized by a pure luminos- 
ity evolution model (PLE) with approximately oc (1 -I- z)^ 
up to z « 2, and consistent with no evolution beyond that 
point. Using a larger ROSAT sample. Page et al. (1997) 
found that PLE underpredicts the number of high-redhsift 
low luminosity AGNs for qo = 0.5. Simple extrapolations 
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of any of the PLE expressions only explain ~ 30 — 50% of 
the soft X-ray Background (0.5-2 keV) by AGN. 

Because of the relatively large PSF of the ROSAT 
PSPC, the identifications of the deepest ROSAT PSPC 
surveys are sometimes ambiguous and misidentifications 
can occur. Based on results of the optical foUowup studies 
of ROSAT VSPG surveys, a number of groups, including 
the Deep ROSAT Smvey (DRS; Griffiths et al. p^96| ) and 



UK Deep Survey (UKD; McHardy et al. |1998D report a 
population of X-ray sources called "Narrow Emission- line 
Galaxies" (NELG) at faint fluxes. On the other hand, faint 
X-ray sources found in the ROSAT Deep Survey on the 
Lockman Hole (RDS-LH), which have accurate source po- 
sitions from 1 million seconds of ROSAT HRI data, are 
still predominantly AGNs down to the faintest fl uxe i n 
the survey (Schmidt et al. 1998; Hasinger et al. 1999| ). 



Some of these have optical spectra which apparently show 
only narrow-lines but have other signs of an AGN activ- 
ity and might have been classified as "NELGs" at the 
criteria of other groups. On the other hand, Lehmann et 
al. ( 1999b| ) have compared redshift distributions of the 
RDS X-ray AGNs, UKD X-ray sources, non X-ray emit- 
ting (at the RDS-LH limit) field galaxies showing narrow 
emission-lines. They found that the redshift distribution of 
UKD X-ray sources has a significant excess over that of the 
RDS-LH sources at z < 0.4. This excess was dominated 
by "NELGs", whose redshift distribution was similar to 
that of non X-ray source narrow emission- line field galax- 
ies. This shows that a significant fraction of "NELGs" 
are likely to be misidentifications by chance coincidences. 
This observation seems to contradict with estimations of 
the relatively low probabilities of such chance coincidences 
by the DRS and UKD groups. A more detailed compari- 
son is urgently needed. Misidentifications affect SXLF es- 
timates in two ways, i.e., by putting a wrong object into 
the sample and by missing the true identifications. Thus it 
is important to have a high spatial resolution image to ob- 
tain unambiguous identifications, especially in the faintest 
regime. 

In this study, we investigate the global behavior of 
the soft X-ray luminosity function (SXLF) of AGNs from 
a combined sample of various ROSAT surveys. We use 
the term "AGN" for both Seyfert galaxies, including type 
I's and type 2's, and QSOs. Preliminary work, using ear- 
lier versions of the combined sample, have been reported 



in Hasinger (1998) and Miyaji et al. (1999a) (hereafter 
M99a), while in this work, we have made a more exten- 
sive analysis with updated ROSAT 'Qvi^i Survey (RBS) 
and ROSAT HccY) Survey (RDS) catalogs including new 
identifcations from observations made in the winter-spring 
season of 1999. In this paper, we put emphasis on the ex- 
pressions representing the global behavior of the SXLF. 
Presenting separate expressions in several redshift inter- 
vals, giving more accurate representation of the data in 
the redshift ranges of interest will be a topic of a future 
paper (Miyaji et al. in preparation, paper II). In paper II, 



we will also present tables of full numerical values of the 
binned SXLF. 

We use a Hubble constant Hq — 50/i5o 
[kms~^ Mpc~^]. The /i5o dependences are explicitly 
stated. We calculate the results with common sets of cos- 
mological parameters: (fimj^A) =(1.0,0.0) and (0.3,0.0). 
For some important parameterized expressions, we also 
show the results for (riin,riA) = (0.3,0.7). 



2. The ROSAT Surveys used in the analysis 

We have used soft X-ray sources identified with AGNs 
with redshift information from a combination of ROSAT 
surveys in various depths/areas from a number of already 
published and unpublished sources. In order to avoid the 
possible bias from the large-scale overdensity and the dis- 
tortion of the redshift-distance relation based on bulk- 



flows in the nearby universe (e.g. TuUy & Shaya 1984), 
we have excluded objects within z < 0.015 from the anal- 
ysis. 

The surveys we have used are summarized in Table nl 
Two optical foUowup programs from the i?OS'^T All-Sky 
Survey (RASS) (Voges 1994), a serendipitous survey of the 
ROSAT PSFC pointed observations (RIXOS), and a num- 
ber of deep pointings specifically aimed for deep surveys. 
Here we describe the AGN sample from each survey. 

All surveys, except for a part of the the Lockman Hole, 
are based on the ROSAT PSPC count rates in the pulse- 
invariant (PI) channel range corresponding to 0.5 - 2 keV. 
For most sources in the Lockman Hole, we have used the 
deeper HRI count rates (see below) with no spectral res- 
olution and sensitive to the 0.1 - 2 keV. 

In order to convert the countrate to flux, we have to 



assume a spectrum. Hasinger et al. ( 1993 ) obtained the 
value of r = 1.96 ± .11 for the average spectral photon in- 
dex i n the L ockman Hole. Other works (Romero-Colenero 
et al. 199C; Almaini et al. 1996| ) also found similar spec- 
tral index for AGNs, but a harder index P « 1.5 for the 
"NELGs". The same class as a part of the population 
they have classified as NELGs may fall into our sample. In 
any case, the ROSAT countrate to the unabsorbed 0.5-2 
keV flux conversion has been made assuming a power-law 
with a photon index of P = 2.0 and corrected for the 
Galactic absorption. In effect, the Galactic column den- 
sity changes the response curves for the fiux-to-countrate 
conversion. However, the extragalactic surveys are mainly 
concentrated on the part of the sky where the Galactic ab- 
sorption is low. Typical values are (0.5 — 1) x 10^*^ [cm~^] 
for the deep surveys and a maximum of 16 x 10^" [cm~^] 
for a small portion of the sky covered by the RBS. Within 
this range, the conversion between the S'x (here and here- 
after, Sx represents the 0.5-2 keV flux and 5x14 is the 
same quantity measured in units of 10^^'*[ergs^^ cm^^]) 
the ROSAT PSPC countrate (in the corresponding chan- 
nel range) only weakly dependent on the spectral shape 
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Table 1. i? 05'^ T" Surveys used in the Analysis 



Survey'' 



[ergs" 



cilim 
"3x14 



Area 

[deg^] 



No. oP 

AGNs 



RBS 


Ri250 


2.0 X 10* 


216 


SA-N 


« 13 


685. 


130 


RIXOS 


3.0 


15. 


205 


NEP 


1.0 


0.21 


13 


UKD 


0.5 


0.16 


29 


RDS-Marano 


0.5 


0.20 


30 


RDS-LH 


0.17-0.9 


0.30 


68 



'^ Abbreviations - RBS: The 7? 05^ T Bright Survey, SA-N: The 
Selected Area-North, RIXOS: The ROSAT International X- 
ray Optical Survey, NEP: The North Ecliptic Pole UKD: The 
UK Deep Survey, RDS-Marano: The ROSAT Deep Survey - 
Marano field, RDS-LH: The ROSAT Deep Survey - Lockman 
Hole. See text for references. ^ Excluding AGNs with z < 0.015. 

and varies by about ±3% for spectral indices T = 2.0±0.7. 
We discuss the conversion for the HRI case in Sect. 2.7. 



For the computation of the SXLF, it is important to 
define the available survey area as a function of limiting 
flux. In case there is incompleteness in the spectroscopic 
identifications, we have made the usual assumption that 
the redshift/classification distribution of these unidenti- 
fied sources is the same as the identified sources at similar 
fluxes. This can be attained by defining the 'effective' sur- 
vey area as the geometrical survey area multiplied by the 
completeness of the identifications. This assumption is not 
correct when the source is unidentified due to non-random 
causes, e.g., no prominent emission lines in the observed 
spectrum. However, given the high completeness of the 
samples used in our analysis, this does not affect the re- 
sults significantly, except for the faintest end of RDS-LH. 
We discuss the effects of the incompleteness at this faint 



end in Sect. 3.5 



Below we summarize our sample selection and com- 
pleteness for each survey. 

2.1. The ROSAT Bright Survey (RBS) 

The RBS program aims for a complete identification of 
the - 2000 brightest sources in the ROSA T All -Sky Sur- 
vey (RASS) for |6| > 30deg (Fischer et al. [1998^ Schwope 
et al. in preparation ) measured in the entire i?05^T band 
(0.1-2.4 keV). For our purposes, we have extracted a sub- 
set selected by the ROSAT Hard band (0.5-2 keV) coun- 
trate of Ci?hard > 0.2 [cts s^^], which makes a complete 
hard countrate-limited sample. Five sources in this sub- 
sample have further been identified as AGNs since M99a 
and included in the analysis. This subsample now has been 
completely identified. 

Since the absorption in our galaxy varies from place 
to place, the same countrate limit corresponds to different 
0.5-2 keV fiux limits based on different galactic A^h values. 
The Nn value range from (0.5 — 16) x 10^° [cm~^]. 
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Fig. 1. The survey area of the combined sample are plot- 
ted as a function of the limiting 0.5-2 keV flux limit (solid 
line). For reference, [A^(> S)]^^ for all the X-ray sources 
is overplotted (dashed-line). 



2.2. The RASS Selected- Area Survey - North 

This survey uses several high galactic latitude areas of 
RASS (a total of 685 [deg^]) for optical identification of 
the sources down to about an order of magnitude fainter 
than the RBS. The fields selected for the survey have the 
Galactic column ranging Nn — (2 — 11) x 10^" [cm^^]. 
Details of the survey have been described in Zickgraf et 
al. ( 1997] ) and the catalog of sourc e iden tifications has 
been published by Appenzeller et al. ( 1998 ). We have fur- 
ther selected our sample such that each field has a com- 
plete ROSAT hard-band (0.5-2 [keV]) countrate-hmited 
sample with complete identifications (Ci?hard > 0.01-0.05 
[ctss-^]). 



2.3. The RIXOS Survey 

The i?OS'^TIntcrnational X-ray Optical Surve y (RI XOS), 
Mason et al. (1999) (see also Page et al. 1996D is a 
serendipitous survey of « 80 PSPC fields covering 15 deg^ 
of the sky. The fiux limit of the deepest field is 5xi4 = 3.0, 
while the actual completeness limit varies from field to 
field. The identification is 97% complete, thus the effect 
of the identification incompleteness is negligible consider- 
ing statistical errors. 

2.4. The North Ecliptic Pole (NEP) Survey 

The data are from Bower et al. ( 1996| ), which gave a cata- 
log of 20 sources in the 15f5 radius region with 5x14 > 1-0. 
One object, RX J1802. 1-1-6629 did not have a redshift en- 
try in Table 2 of Bower et al., but in the text, they argued 
that the most probable interpretation of this object was 
a weak- lined QSO at z ^^ 1. Thus we have assigned a 
redshift of 1.0 to this source. There is one unidentified 
source, making the idcntificaton of the sample 95% com- 
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plete. Thus we have set the effective survey area of the 
NEP survey as 95% of the geometrical area. 

2. 5. The UK Deep Survey 

Based on a 115 [ks] of ROSAT PSPC observation, 
McHardy et al. (1998) pubhshed a hst of sources and iden- 
tification of X-ray sources down to 5x14 = 0.19. A signif- 
icant fraction of their identifications are "NELGs" (Nar- 
row Emission-Line Galaxies) and the fraction increases 
towards fainter fluxes. As mentioned in Sect. H^, a part of 
these NELGs are likely to be misidentifications. The iden- 
tifications of other NELGs might be correct, but those 
would have been classified as AGNs with the criteria of 
Schmidt et al. (|1998D. 



To include the results of the UKD survey in our sam- 
ple, we would like to include their NELGs in the latter 
category, but exclude those in the former category. We 
find that the redshift distribution of the AGN-hNELG 
classes in the UKD survey is significantly different from 
the AGN-|-galaxy classes in the Lockman survey if we in- 
clude all sources down to 5x14 > 0.19. If we limit the 
sample to brighter sources (5xi4 > 0.5), the redshift dis- 
tributions are consistent with each other. Thus, in this 
work, we limit the samples from UKD and other deep 
PSPC surveys to 5xi4 > 0.5, assuming that the misidenti- 
fication problem would not affect the analysis significantly 
above this limit. 

2.6. The ROSAT Deep Survey - Marano Field 

For the same reason as the UKD case, we have also used 
the same flux cutoff 5xi4 > 0.5 for the survey in the 15'- 
radius region on the Marano field (Zamorani et al. 1999| ), 
based on a deep PSPC exposure. Source fluxes of their cat- 
alog have been updated since the version used by M99a. 
The identiflcations are 100% complete for the 14 sources 
in 5x14 ^ 1 and 4 of the 27 sources remain unidentifled 
or ambiguous (85% complete) in 0.5 < 5xi4 < 1. As be- 
fore, we have reduced the survey area by 15% in this flux 
range to deflne the effective survey area used in the SXLF 
calculations. 



2.7. The ROSAT Deep Survey - Lockman Hole 

There are 200 ks of PSPC and 1 Msec of HRI observations 
on this field (Hasinger et al. 



1998) 



The source list and 
identifications for the brightest 50 sources (5xi4 > 0.5) 
have been published (Schmidt ct al. 1998[ ). In this work, 
we have included further unpublished identifications down 
to 5x14 — 0.17. These include identifications and redshifts 
based on spectra obtained in March 1998 with the Keck 



10m telescope (Hasinger et al. 1999), which have also been 



The conversion between the HRI countrate and the 
0.5-2 keV flux has been determined from the mean values 
of overlapping sources between the HRI and PSPC. The 
convsersion carries more uncertainties based on spectra, 
because the HRI has practically no spectral resolution and 
has some sensitivity down to 0.1 keV. With the HRI, the 
conversion factor varies by ±40% for photon indices P = 
2.0 ±0.7. 

The basic strategy of defining the combined PSPC- 
HRI sample has been explained in Hasinger et al. (1999). 
In this paper, we have slightly modified the fiux-limit and 
areas of the HRI sample in order to optimize our AGN 
sample in the presence of new identifications: 

— We use the deeper HRI-detected sample and HRI 
fluxes for the region 12.0 arcminutes from the HRI cen- 
ter (0.126 deg2). At the faintest fluxes (0.17 < 5xi4 < 
0.24), we have further limited the area to 10.1 arcmin- 
utes from the HRI center (0.090 deg^). This choice al- 
lows us to avoid the problem of incomplete source de- 
tection due to source confusion (see Fig. 2b of Hasinger 
et al. |1999D . A total of 48 AGNs are present in this HRI 
sample. 

— Outside of the HRI region dcflned above, and within 
18.4 arcminutes from the PSPC center, we have used 
the PSPC detected sources and PSPC fluxes. This cor- 
responds to 0.175 deg^. For completeness, we have im- 
posed a flux cutoff of 5x14 > 0.38 for PSPC off-axis 
angles smaller than 12'.5 and 5xi4 > 0.97 for PSPC 
offaxis angles between 12.5 and 18.4 arcminutes re- 
spectively. 

— The sources in the HRI/PSPC combined sample have 
been 100% identifled for 5xi4 > 0.38. Four of the 31 
sources in 0.17 < 5xi4 < 0.38 remain spectroscopically 
unidentifled. Thus we have reduced the effective survey 
area by 13% for 0.17 < 5xi4 < 0.38 to compensate for 
the identiflcation incompleteness. 



2.8. The combined sample 

In our combined sample, there are 691 AGNs ranging from 
4.2x10"" to 1.7xl0~i^[ergs-^cm-2]. The effective sur- 
vey area for the combined sample is plotted as a func- 
tion of the limiting flux in Fig. ^, overlaid with the value 
of A^(> 5)~^, showing that the combined sample indeed 
covers the above flux range continuously. 

The redshift - luminosity scatter digrams of the sample 
objects are shown in Fig. for the (fJoi^A) = (1.0,0.0) 
universe. In any case, the luminosities have been calcu- 
lated by: 



i. 



4^dL(^)'5o.5- 



2kcV 



(1) 



included in M99a. Further four spectroscopic identiflca- 
tions obtained with the Keck telescope in February 1999 
have been added to the sample since M99a. 



where di^{z) is the luminosity distance as a function of 
redshift, which depends on the choice of cosmological pa- 
rameters. This corresponds to the no K-correction case. 
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Fig. 2. The AGNs in the combined sample are plotted in 
the z — LogL^ for (SljilA) =(1.0,0.0). Different symbols 
correspond to different surveys as labeled. 

Explanation on our K-correction policy is explained in 



Sect. 3.1. Hereafter, the symbol L^ refers to the quan- 
tity defined in Eq. (|^) expressed in units of /i^g^[ergs~^], 
unless otherwise noted. 



3. The ROSAT AGn SXLF 

3.1. K-Correction and AGN subclasses 

In this section, we choose to present the SXLF in the ob- 
served 0.5-2 keV band, i.e., in the 0.5(1 + z)- 2(1 + z) keV 
range in the object's rest frame. Thus no K-correction 
has been appplied for our expressions presented in this 
section. Also we choose to include all emission-line AGNs 
(i.e., except BL-Lacs), including type I's and type 2's. The 
primary reason for these choice is to separate the model- 
independent quantities, directly derived from i?OS'^T sur- 
veys, from model-dependent assumptions. Here we explain 
the philosophy behind these choices in detail. 

There are a variety of AGN spectra in the X-ray 
regime, but the information on exact content of AGNs 
in various spectral classes is very limited. Currently pop- 
ular models explaining the origin of the 1-100 keV CXRB 
involve large contribution of self-absorbed AGNs (Madau 
et al. |l994|; C omastri et al. |1995| Miyaji et al. |1999b| ; Gilli 
et al. 1999). Although they are selected against in the 
ROSAT band, some of these absorbed AGNs come into 



our sample. These absorbed AGNs certainly have different 
K-correction properties than the unabsorbed ones. While 
these absorbed AGNs are mainly associated with those 
optically classified as type 2 AGNs, the correspondence 



between the optical classification and the X-ray absorp- 
tion is not straightforward. Especially, there are many 
optically type-1 AGNs (with broad-permitted emission 
lines), which show apparent X-ray absorption of some 
kind. For example, a number of Broad Absorption Line 
(BAL) QSOs are known to have strongly absorbed X-ray 
spectra (e.g. Mathur et al. 1995; Gallagher et al. 1999). 



At the fainter/high-redshift end of our survey, there may 
be some broad-line QSOs of this kind or some intermedi- 
ate class. Broad-line AGNs with hard X-ray spectra have 
been found in a number of hard surveys (Fioreet al. 1999; 
Akiyama et al. |l999| ). In Schartel et al. (|l997|) 's study, 
all except two of the 29 AGNs from the Piccinotti et al.'s 



( |1982D catalog have been classified as type I's, but about a 
half of them show X-ray absorption, some of which might 
be caused by warm absorbers. In view of these, using only 
optically-type 1 AGNs to exclude self-absorbed AGNs is 
not appropriate. Also optical classification of type 1 and 
type 2 AGNs depend strongly on quality of optical spec- 
tra. Thus classification may be biased, e.g. as a function 
of flux. However, the SXLF for the type 1 AGNs is of his- 
torical interest and shown in Appendix A. As shown in 
Appendix A., non type-1 AGNs are very small fraction of 
the total sample and excluding these does not change the 
main results significantly. 

On the other hand, our sample of 691 AGNs with 
extremely high degree of completeness carries little un- 
certainties in the fluxes in the 0.5-2 keV band in the 
observer's frame, redshifts, and classiflcation as AGNs. 
Thus, we choose to show the SXLF expression in the ob- 
served 0.5-2 keV band, or 0.5(l-|-z)-2(l-|-z) keV band at 
the source rest frame, in order to take full advantage of this 
excellent-quality sample without involving major sources 
of uncertainties. The expressions in the observed band 
may have less direct relevance for discussion on the ac- 
tual AGN SXLF evolution. However they are more useful 
for discussing the contribution of AGNs to the Soft X-ray 
Background (Sect. ||), interpretation of the fluctuation of 
the soft CXRB, and evaluating the selection function for 
studying clustering properties of soft X-ray selected sam- 
ple AGNs. 

In practice, the expressions can also be considered a 
K-corrected SXLF at the zero-th approximation, since ap- 
plying no K-correction is equivalent to a K-correction as- 
suming r = 2. This index has been historically used in pre- 



vious works (e.g. Maccacaro et al. 1991; Jones et al. 1996), 
thus our expression is useful for comparisons with previous 
results. A r = 2 power-law spectrum can be considered 
the best-bet single spectrum characterizing the sample, 
because in the i? 05^1" sample, absorbed AGNs (includ- 
ing type 2 AGNs, type 1 Seyferts with warm absorbers, 
BAL QSOs) are highly selected against. Nearby type 1 
AGNs show an underlying power-law index of F = 2 at 
E Z l[keV] (e.g. George et al. |1998D , which is the en- 
ergy range corresponding to 0.5-2 keV for the high red- 
shifts where K-correction becomes important. The reflec- 



Miyaji, Hasinger, Schmidt: ROSAT AGf^ Luminosity Function 



tion component, which makes the spectrum apparently 
harder, becomes important only above 10 keV. This is 
outside of the ROS AT hand even at z > 4. The above ar- 
gument is consistent with the fact that the average spectra 
of the faintest X-ray sources, especially those indentified 
with broad- line AGNs, have F sa 2 (Hasinger et al. 1993; 



Romero- Colmenero et al. 1996; Almaini et al. 1996) in 



the ROS AT band. Therefore, at the zero-th approxima- 
tion, one can view our expression as a K-corrected SXLF 
of AGNs, especially at high luminosities. The goodness of 
this approximation is highly model-dependent and a dis- 
cussion on further modeling beyond this zero-th approxi- 
mation is given in Sect. |6|. 

3.2. The bmned SXLF of AGNs 

The SXLF is the number density of soft X-ray-selected 
AGNs per unit comoving volume per Log L^ as a function 
of Lx and z. We write the SXLF as: 



d$ 



d Log Lx 



(Log L^,z). 



Fig.H shows the binned SXLF in different redshift 
shells estimated using the ^ 1/14 estimator: 



d$ 



d Log Lx 



( Log Lx 



(A Log Lx)j 



(2) 



where the L^ — z bins are indexed by j and AGNs 
in the sample falling into the j'-th bin are indexed by i, 
Vi{L^) is the available comoving volume in the redshift 
range of the j — th bin where an AGN with luminosity 
Lx would be in the sample. The luminosity function is 
estimated at (Log Lxj,'Zj), where a bar represents the V'-^ 
weighted average over the AGNs falling into the j-th bin. 
Also (A Log Lx)j is the size of the j— th bin in Log L^. 

Rough Ict errors have been estimated by: 



d$ 



d Log Lx 



(Log Lx 



(A Log Lx)j 



(3) 



In case there is only one AGN in the bin, we have plotted 
error bars which correspond to the exact Poisson errors 
corresponding to the confidence range of Gaussian la. In 
this way, we can also avoid infinitely extending error bars 
in the logarithmic plot. 

Fig. ^(a)(b) shows the binned SXLF calculated for 
(17„i,fiA) = (1.0,0.0) and (0.3,0.0) respectively. In Fig. |, 
we have also plotted some interesting upper- limits, in case 
there is no object in the bin. In the figure, we show upper 
limits corresponding to 2.3 objects (90% upper-limit). See 
caption for details. 

We note that the binned ^ 1/14 estimate can cause 
a significant bias, especially because the size of the bins 
tend to be large. For example, at low luminosity bins with 
corresponding fluxes close to the survey limit, the value 
of T4 can vary by a large factor within one bin. Also the 



choice of the point in L^ space representative of the bin, 
at which the SXLF values are plotted, may change the im- 
pression of the plot significantly. Thus the SXLF estimates 
based on the binned J2 1/^a can be used to obtain a rough 
overview of the behavior, but should not be used for sta- 
tistical tests or a comparison with models. Full numerical 
values of the binned SXLF including X) 1/^ values, im- 
proved estimations by a method similar to that discussed 



by Page & Carrera 1999, and the numbers of AGNs in 



each bin will be presented in paper II. 

A number of features can be seen in the SXLF. As 
found previously, our SXLF at low z is not consistent with 
a single power-law, but turns over at around Log L^ ^ 
43 — 44. The SXLF drops rapidly with luminosity beyond 
the break. We see a strong evolution of the SXLF up to the 
0.8 < z < 1.6 bin, but the SXLF does not seem to show 
significant evolution between the two highest redshift bins. 
Figs. H (a) (b) show that these basic tendencies hold for the 
two extreme sets of cosmological parameters. 

3.3. Analytical expression - statistical method 

It is often convenient to express the SXLF and its evolu- 
tion in terms of a simple analytical formula, in particular, 
when using as basic starting point of further theoretical 
models. 

Here we explain the statistical methods of parameter 
estimations and evaluating the acceptance of the models. 
A minimum x^ fitting to the binned 1/T4 estimate is not 
appropriate in this case because it can only be applied to 
binned datasets with Gaussian errors and at least 20-25 
objects per bin are required to achieve this. In our case, 
such a bin is typically as large as a factor of 10 in L^ 
and a factor of two in z, thus the results would change 
depending where in the [z Lx) bin the comparison model 
is evaluated. 

The Maximum-Likelihood method, where we exploit 
the full information from each object without binning, 
is a useful method for parameter estimations (e.g. Mar- 
shall et al. 198S), while, unlike x^, it does not give abso- 
lute goodness of fit. The absolute goodness of fit can be 
evaluated using the one-dimensional and two-dimensional 
Kolgomorov-Smirnov tests (hereafter, ID-KS and 2D-KS 
tests respectively; Press et al.l992; Fasano & Franceschini 



1987) to the best-fit models. 



As our maximum-likelihood estimator, we define 



C 



-2 El- 



iV(LogLxj,Zi) 



/ / N{Log Lx, z)d LogLx dz 



(4) 



where i goes through each AGN in the sample and 
N{Log Lx, z) is the expected number density of AGNs 
in the sample per logarithmic luminosity per redshift, cal- 
culated from a parameterized analytic model of the SXLF: 



Miyaji, Hasinger, Schmidt: 7?OS'j4r AGN Luminosity Function 






com 



O 



t3 



10-4 
10-5 

10-6 
10-^ 

10-8 
10-9 

10-1° 
10-11 

10-12 



I I I I I I I I 



(a) (n^,nj=(i. 0,0.0) 




O z:0.4-0 
■ z:0.2-0.4 
□ z:0.015-0.2 



I I I I I I 



I I I I I I I I 



± 



± 



(b) (n ,nJ = (o.3,o.o) 




42 43 44 45 46 47 
Log L^ [h^s erg s-i] 



42 43 44 45 46 47 
Log L^ [hg^ erg g-i] 



Fig. 3. The ^V~^ estimates of the SXLFs are plotted with estimated la errors. Different symbols correspond to 
different redshift bins as indicated in the panel (a) and data points belonging to the same redshift bin are connected. 
The position of the symbol attached to a downward arrow indicates the 90% upper limit (corresponding to 2.3 objects), 
where there is no AGN detected in the bin. 



d Log L, 



N(hogL^,z) = 
■dA{z)^{l + zf 



Except for the global normalization A, we have made 



• dz 



{z)-A{L^/dl), 



(5) 



where dA{z) is the angular distance, dT/dz[z) is the dif- 
ferential look back time per unit z (e.g. Boldt 1987) and 
A[S-x) is the survey area as a function of limiting X-ray 
flux (Fig. El). Minimizing C with respect to model param- 
eters gives the best-fit model. Since A£ from the best-fit 
point varies as Ax^, we determine the 90% errors of the 
model parameters corresponding to A£ = 2.7. The min- 
imizations have been made using the MINUIT Package 
from the CERN Program Library (James 1994). 

Since the likelihood function Eq. (0) used normalized 
number density, the normalization of the model cannot be 
determined from minimizing C, but must be determined 
independently. We have determined the model normaliza- 
tion (expressed by a parameter A in the next subsections) 
such that the total number of expected objects (the de- 
nominator of the right-hand side of Eq. (^) is equal to 
the number of AGNs in the sample {N°^''^). 



use of the MINUIT command MINOS (see James |1994D 
to serach for errors. The command searches for the pa- 
rameter range corresponding to A>C < 2.7, where all other 
free parameters have been re-fitted to minimize C dur- 
ing the search. The estimated 90% confidence error for A 
is taken to be \.7A {\/W°^)~^ and does not include the 
correlations of errors with other parameters. 



The ID-KS tests have been applied to the sample dis- 
tributions on the Lx and z space respectively. The 2D- 
KS test has been made to the function N{Log L^,z). 
We have shown the probability that the fitted model is 
correct based on the ID- and 2D-KS tests. For the 2D- 
KS test, calculated probability corresponding to the D 
value from the analytical formula is accurate when there 
are ^ 20 objects and the probabilities ^ 0.2. If we ob- 
tain a probability ^ 0.2, the exact value does not have 
much meaning but implies that the model and data are 
not significantly different and we can consider the model 
acceptable. We have searched for models which have ac- 
ceptance probabilities greater than 20% in all of the KS 
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tests. Strictly speaking, the analytical probability from 
the KS-test D values are only correct for models given a 
priori. If we use paramters fitted to the data, this would 
overestimate the confidence level. A full treatment should 
be made with large Monte-Carlo simulations (Wisotzki 



Table 2. Best-fit PLE and PDE Parameters 



199S), where each simulated sample is re-fitted and the 



D- value is calculated. However, making such large simula- 
tions just to obtain formally-correct probability of good- 
ness of fit is not worth the required computational task. 
Instead, we choose to use the analytical probability and 
set rather strict acceptance criteria. 

3.4- Analytical expression - overall AGN SXLF 

Using the method described above, we have searched for 
an analytical expression of the overall SXLF. The overall 
fit has been made for the redshift range 0.015 < z < 5. 
Also for the fits, we have limited the luminosity range to 
Log L^ ^ 41.7. 

As described in Sect. ||, the lower redshift cutoff is 
imposed to avoid effects of local large scale structures, 
which may cause a deviation from the mean density of 
the present epoch and thus can cause significant bias to 
the low luminosity behavior of the SXLF. At the lowest 
luminosities (Log L^ < 41.7), there is a significant excess 
of the SXLF from the extrapolation from higher lumi- 
nosities. This excess conn ects w ell with the nearby galaxy 
SX LF by Schmidt et al. ( |l996| ) (see also e.g. Hasinger et 
al. 1999) and may well contain contamination from star 



formation activity (see also Lehmann et al. 1999a ). For 
finding an analytical overall expression, we have not in- 
cluded the AGNs belonging to this regime. 

As an analytical expression of the present-day (z — 
0) SXLF, we use the smoothly-connected two power-law 
form: 

^ "d^Log Z "^ = ^ [(^'^z^*)" + iL./L..rr' (6) 

As a description of evolution laws, the following models 
have been considered: 



3.4.1. Pure-luminosity and pure-density evolutions 



As so me pr evious works (e.g. D elia Ceca et al. |1992| ; Boyle 
et al. 1994 : Jones et al. 1996 ; Page et al. 1996 ), we have 
first tried to fit the SXLF with a pure-luminosity evolution 
(PLE) model. 



d$(Xx,2) _d$(L,/e(z),0) 



d Log Lx 



d Log Lx 



(7) 



For the evolution factor, wc have used a power-law 
form: 



eiz) 



(l + z)Pi iz<z,) 

e(z,)[(l + z)/(l + Ze)F2 (z>z,) 



(8) 



Model 

{flm, ^a) 



Parameters^ /KS probabilities 



PLE 
(L0,0.0) 



PLE 
(0.3,0.0) 



PDE 

(L0,0.0) 



PDE 

(0.3,0.0) 



A = (4.0 ± .3) X 10""; L, = 0.33 ± .10 

71 =0.60 ±.16; 72 = 2.34 ± .12;pl = 3.0 ± .2 

Zc = 1.42 ±. 17; p2 = 0.31?;^ 

Pks = .002, 3 10~^ 1 10-^ (for L,z,2D) 



A = (3.1 ± .2) X 10"^ L, = 0.38 ± .12 

71 =0.57 ±.16; 72 = 2.35±.12;pl = 2.9 ± .2 

Zc = 1.54±.25;p2 = 0.3±.7 

Pks = .08, .001,2 10"^ (for L,z,2D) 

A = (6.0 ± .4) X 10"^ L, = 1.08 ± .4 

71 =0.74 ±.13; 72 = 2.28±.ll;pl = 4.6 ± .3 

Zc = 1.60±.25;p2 = 0.6±l.l 

Pks = 0.9,0.9,0.16 (for L,z,2D) 

A = (5.4 ± .3) X lO"'^; L, = 1.13 ± 0.4 

71 =0.76 ±.13; 72 = 2.22 ± .10; pi = 4.6 ± .3 

Zc = 1.62±.26;p2 = 1.3 ± 1.1 

Pks = 0.8,0.7,0.1 (for L,z,2D) 



^'Units - A: [hlo Mpc-»], L,: [lO"** /igo'erg s'^], /xi2: 
[10~ ergs~^cni~ deg~ ] in 0.5-2 keV. Para mete r errors cor- 
respond to the 90% confidence level (see Sect. 3.3). 



The best-fit values are listed in the upper part of Table 
H along with ID-KS and 2D-KS probabilities using the 
analytical formula. In Table g and later tables, the three 
values of Pks represent the probabilities that the model is 
acceptable for the ID-KS test in the L^ distribution, ID- 
KS test in the z distribution, and 2D-KS test in the {L^,z) 
distribution respectively. Note that there are cases which 
are accepted by ID-KS tests in both distributions but fail 
in the 2D-KS test. The results of the fit show that the 
PLE model is certainly rejected with a 2D-KS probability 
of P2DKS = 5 X IQ-'^ and 1 x 10"^ for the ^^^=1 and 0.3 
(flA = 0) cosmologies respectively. 

As an alternative, we have also tried the Pure-Density 
Evolution model (PDE), which seemed to fit well in our 
preliminar y an alysis for the fl^—1 (fljy — 0) universe 
(Hasinger ^IH). 



d$(Lx,z) ^ d^(i:,/,0) 
d Log Lx d Log Lx 



eiz) 



(9) 



where e{z) has the same form as Eq. m). The 2D-KS 
probabilities are P2DKS =0.16 and 0.1 for the flm—^ and 
0.3 (flA — 0) respectively. Thus the acceptance of the over- 
all fit is marginal, especially for riin=l- However the PDE 
model has a serious problem of overproducing the soft X- 
ray background (Sect. ^. For a further check, we have 
made separate fits to high luminosity (Log L^ > 44.0) 
and low luminosity (Log Lx > 44.0) samples to compare 
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Table 3. Best-Fit LDDEl Parameters 



Model 



Parameters^/KS probabilities 



LDDEl 
(1.0,0.0) 



LDDEl 

(0.3,0.0) 



LDDEl 

(0.3,0.7) 



A = (1.01 ± .06) X 10"''; L. = 0.75^*;; 

7i =0.75 ±.15; 72 = 2.25±.10;pl = 5.1 ± , 

2c = 1.57 ± .15; p2 = 0.0 (fixed) 

a = 1.7 ± .8; LogLa = 44.1 (fixed) 

Pks = 0.6,0.4,0.5 (for L,2,2D); 



;L*=0.56t;?g 



A^ (1.56 ±.10) X 10" 

71 =0.68 ±.18; 72 = 2.19 ±.08; pi = 5.3 ± .4 

2c = 1.59 ± .14; p2 = 0.0 (fixed) 

a = 2.3 ± .7; LogLa = 44.3 (fixed) 

Pks = 0.5,0.3,0.3 (for L,2,2D) 

A = (1.61 ± .10) X 10-"; L, = 0.56tg 

71 =0.66 ±.18; 72 = 2.19 ±.08; pi = 5.3 ± .4 

2c = 1.58 ± .14; p2 = 0.0 (fixed) 

a = 2.6 ± .7; LogLa = 44.4 (fixed) 

Lks = 0.4,0.4,0.3 (for L,2,2D) 



^Units - A: [/lio Mpc"^], L*: [lO''"' /i^o^erg s % Parameter 



errors correspond to the 90% confidence level (see Sect. 3.3) 



the evolution index pi in 0.015 < z < 1.6 for Vt^^ = 1.0. 
We have obtained pi = 5.3 ± 0.5 and 4.1 ± 0.5 (90% er- 
rors) for the high and low luminosity samples respectively. 
Thus the density evolution rate is somewhat slower at low 
luminosities. Of course at the low luminosity regime, the 
fit was weighted towards nearby objects. If the evolution 
does not exactly follow the power-law form (ex [1 ± z]p^), 
spurious difference in evolution rate can arise. Visual in- 
spection of Fig. ^ might suggest that at z < 0.4, the evo- 
lution rate seems larger at low luminosities, as opposed to 
the results shown above for z < 1.6. However, perform- 
ing the same experiment for the z < 0.4 AGNs showed 
pi = 5.7± 1.8 and 5.8± 1.2 for the high and low luminos- 
ity samples respectively, indicating no difference within 
relatively large errors. For the 0.4 > z > 1.6 sample, the 
results are pi — 6.2 ± 0.8 and 3.0 ± 1.0, again, for the 
high and low luminosity samples respectively. This differ- 
ence and the soft CXRB overproduction problem lead us 
to explore a more sophisticated form of the overall SXLF 
expression as described in the next section. 

3.4.2. Luminosity-dependent density evolution 

We have tried a more complicated description by modify- 
ing the PDF model such that the evolution rate depends 
on luminosity (the Luminosity-Dependent Density Fvolu- 
tion model). In particular, as shown above, it seems that 
lower evolution rate at low luminosities than the PDF 
case would fit the data well. This tendency is also seen 
in the optical luminosity function of QSOs (Schmidt & 



Green 1983| ; Wisotzki |l998| ). The particular form we have 
first tried (the LDDFl model) replaces e(z) in Eq. (m by 
e{z, Lx), where 



e{z,L^) = 

(1 + 2)'°a'^(0'Pl-"Log [La/Lx]) (2 < z^; Lx < La) 

(l + z)Pi {z<z,;L^>L,) 

(Z > Zc) 



{l + z)P^ 

e(Ze,ix)[(l±z)/(l±^c)]^' 



(10) 



In Fq. (pi)|). The parameter a represents the degree 
of luminosity dependence on the density evolution rate 
for Lx < La. The PDF case is a = and a greater value 
indicates lower density evolution rates at low luminosities. 

The best-fit LDDFl parameters, the results of the KS 
tests, and the integrated 0.5-2 keV intensity are shown in 
Table |[ Table ^ shows that considering the luminosity 
dependence to the density evolution law has significantly 
improved the fit. The 2D-KS probabilities (analytical) are 
more than 30% for all sets of cosmological parameters. 

We have considerd another form of the LDDF model 
(designated as LDDE2), which was made to produce 90% 
of the estimated 0.5-2 keV extragalactic background. The 
details of the construction of the LDDF2 is discussed in 
Sect. Q where the contribution to the Soft Cosmic X-ray 
Background is discussed. In figures in the following dis- 
cussions, the LDDF2 model is also plotted. 

For an illustration, in Fig. ^ we show the behavior of 
the density evolution index for z < Zc as a function of lu- 
minosity for our PDF, LDDEl and LDDF2 models. Fig | 
shows the behavior of the model SXLFs at z=0.1 and 1.2. 
In this figure, only the part drawn in thick lines is con- 
strained by data and thin lines are model extrapolations. 
These figures are only meant for illustrative purposes and 
thus are only shown for the (O^, riA)=(0.3, 0) cosmology, 
where differences among models are more pronounced. 

3. 5. Comparison of the data and the models 

For a demonstration of the comparison between the an- 
alytical expressions and the data, we have plotted the 
5'i-5iV(> S) curve (the Log N - Log S curve plotted 
in such a way that the Euclidean slope becomes horizon- 
tal) for AGNs in our sample with expectations from our 
models (Fig. g). Also the redshift distribution of the sam- 
ple has been compared with the models in Fig. |^. These 
two comparisons already show intersting features. As ex- 
pected, the PLF underpredicts and PDF overpredicts the 
number counts of lowest flux sources. In the redshift dis- 
tribution, the PLE overpredicts the number of z ^ 0.08 
sources while it slightly underpredicts the z « 1 sources. 
Although the deviation in each redshift bin seems small, 
the deviations in the neighboring bins are consistent and 
these systematic deviations can be sensitively picked up 
by the KS test in the z distribution (see small values of 
the Pks in z for the PLE model in Table |). 
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Fig. 4. The behavior of the evolution indices a,t z <i Zc are 
shown as a function of luminosity for various density evo- 



lution models : PD E (short-dashed, Sect. |3.4.l| ), LDDEl 
(long-dashed |X|), and LDDE2 (dot-dashed, Sect. |). 
The lines for the (r2ni,i^A) = (0.3,0) case are shown. 
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Fig. 5. The behavior of the model SXLFs at z=0.1 and 1.2 
are shown respectively for the PLE (dotted) , PDE (short- 
dashed), LDDEl (long-dashed), and LDDE2 (dot-dashed) 
models. For the z=1.2 curves, thick-line parts show the 
portion covered by the sample (5'xi4 ^ 0.2) and the thin- 
line parts are extrapolations to fainter fluxes. The lines 
are for {nm,^A) = (0.3,0). 

The plots in Figs. ^ and |^ are comparisons of distribu- 
tions in one-dimensional projections of a two-dimensional 
distribution. Only with these projected plots, one can eas- 
ily overlook important residuals localized at certain loca- 
tions. Thus we also would like to show the comparison 
in the full two-dimensional space. In literature, models 
are often overplotted to the binned SXLF plot calculated 
by the V~^ estimate like Fig. ^. However, given unavoid- 
able biases associated with the binned V~^ estimates (see 
Sect. 3^), such a plot can cause one to pick up spurious 
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S^ (0.5-2 keV) [erg crri-^ s-i] 

Fig. 6. The S^^^N{> S) (a horizontal line corresponds 
to the Euclidean slope) curve for our sample AGNs is 
plotted with 90% errors at several locations and are com- 
pared with the best-fit PLE (dotted), PDE (shot-dashed), 
LDDEl (long-dashed), LDDE2 (dot-dashed) models for 
the (fimj^A) = (IjO) (upper panel) and (0.3,0) (lower 
panel). The thin-solid fish is from the fluctuation analysis 
of the Lockman Hole HRI data (including non-AGNs) by 
Hasinger et al. (1993) 



residuals. Thus we have plotted residuals in the following 
unbiased manner. For each model, we have calculated the 
expected number of objects falling into each bin (7V'°°'^°') 
and compared with the actual number of AGNs observed 
in the bin (TV^ata-j^ r^^iQ full residuals in term of the ratio 
^data/^modei ^^.g plotted in Fig. | for the PDE, LDDEl 
and LDDE2 models for two sets of cosmological parame- 
ters as labeled. The error bars correspond to la Poisson 
errors (cp) estimated using Eqs. (7) and (11) of Gehrels 
( 1986| ) with 5 = 1. Points belonging to different redshift 
bins are plotted using different symbols as labeled (iden- 
tical to those in Fig. y). These residual plots show which 
part of the z — Lx space the given models are most rep- 
resentative of, which part is less constrained because of 
the poor statistics, and where there are systematic residu- 
als. It seems that the models underpredict the number of 
AGNs in the highest luminosity bin at 2.3 < z < 4.6 by a 
factor of 10, but statistical significance of the excess is still 
poor (2 objects against the models predictions of about 
0.2). These AGNs do not constrain the fit strongly and 
excluding them did not change the results significantly. 
Also there is a scatter up to a factor of 2 from the model 
in 45 Si LogLx SS 46, but no points are more than 2a 
away from either of the LDDEl and LDDE2 models in 
both cosmologies. 

The only data point which is more than 2a away from 
LDDEl or LDDE2 model is the lowest luminosity bin at 
1.6 < z < 2.3 (filled triangle), i.e., 43.6 <, LogLx ^ 44.2 
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Fig. 7. The redshift distribution of the AGN sample, his- 
togrammed in equal interval in log z, is compared with 
predictions from the best-fit PLE (dotted), PDE (shot- 
dashed), LDDEl (long-dashed), and LDDE2 (dot-dashed) 
models for two sets of cosmological parameters as labeled. 
The assymmetric error bars correspond to approximate 
Icr Poisson errors calculated using Eqs. (7) and (11) of 
Gehrels (Il98^) with 5=1. 



for (f7„i,r2A) = (1.0,0.0) or 43.8 <, LogL^ <, 44.5 for 
(fim.^A) = (0.3,0.0). Both LDDEl and LDDE2 models 
overpredict the number of AGNs by a factor of ~ 2 in both 
cosmologies, which are 2.2 — 3.8a away. However, this lo- 
cation corresponds to the faintest end of the deep surveys 
with a certain amount of incompleteness in the identifica- 
tions. Our incompleteness correction method (Sect. ||) is 
valid only if the unidentified source are random selections 
of the X-ray sources in the similar flux range. However, 
these sources have remained unidentified because of the 
difficulty of obtaining good optical spectra and not by a 
random cause. Thus it is possible that the incompleteness 
preferentially affects a certain redshift range. Actually the 
deficiencies were much larger in the previous version (see 
Fig. 8 of Hasinger et al. 1999 ). The discrepancies decreased 
after the February 1999 Keck observations of the faintest 
Lockman Hole sources with rather long exposures, where 
three of the four newly identified source turned out to 
be concentrated in this regime. Thus it is quite possible 
that the remaining four unidentified sources are also con- 
centrated in this regime. In that case, the LDDE mod- 
els can also fit to this bin within 2a. Actually the newly 
identified and unidentified sources typically have very red 
R—K' colors (Hasinger et al. 1999; Lehmann et al. 1999q), 
which probably belong to a similar class to those found by 
Newsam et al. (1998). If the red R — K' color comes from 



the stellar population of underlying galaxy, they are likely 
to be in a concentrated redshift regime. On the other hand, 
if it represents obscured AGN component, they can be in 
a variety of redshift range. At this moment, it is not clear 
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Fig. 8. The full residuals of the fit are shown for the 
PDE, LDDEl and LDDE2 models in two sets of cosmo- 
logical parameters as labeled in each panel. The redid- 
ual in each bin has been calculated from actual number 
of sample AGNs falling into the bin and the model pre- 
dicted number. Different symbols correspond to different 
redshidt bins as indicated above the top panel, which are 
identical to those used in Fig.g. One sigma errors have 
been plotted using approximations to the Poisson errors 
given in Gehrels ( |1986| ). The upper limit corresponds 2.3 
objects (90% upper-limit). 



whether the deficiencies in this location is due to incom- 
pleteness or indicate an actual behavior of the SXLF. 

Based on the results of the 1-D and 2-D KS tests, we 
have rejected the PLE model. We favor the LDDEl and 
LDDE2 models over the PDE model based on the KS 
tests and as well as the CXRB constraints (see below). It 
may be interesting to show the exact location where the 
largest discrepancies are for these models, as compared 
to the LDDE models. This can be most clearly shown 
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Fig. 9. Residuals in the -x space (see text) are shown for 
two resdhift bins, i.e., 0.015< z < 0.2 and 0.8< z < 1.6, 
where differences among different models are apparent. 
Different line styles correspond to different models. See 
caption for Fig. ra for the line styles. The luminosity bins 
are shown as horizontal bars bordered by ticks. 



by plotting residuals in the -x = {N'^°'^°^ - N'^''^'')/a^ 
space. We have shown the — x residuals for redshift bins 
where there are notable differences among these models, 
i.e., 0.015 < z < 0.2 and 0.8 < z < 1.6. These are shown 
in Fig. 0. For both cosmologies, the PLE model systemat- 
ically shifts from overprediction to underprediction with 
increasing luminosity at the lowest redshift bin. At the 
higher resdhift bin, the opposite shift can be seen. The 
curve converges closer to zero at both high and low lumi- 
nosity ends just because there are only small numbers of 
objects in these bins causing poor statistics. More appar- 
ently in the (r2m,riA) — (0.3,0) universe, the PDE model 
also shows a significant scatter around zero. 

The data in the lower luminosity part 42 ^ LogLx ~ 
43.5 in the lowest redshift bin (0.015 < z < 0.2) are 
crucial in rejecting the PLE model, as seen in Figs. ^ 
and y. This regime, consisting of '^ 90 AGNs, has low 
SXLF values compared with the PLE extrapolation from 
the higher redshift data. Actually we cannot discriminate 
between the PLE and LDDE models for the sample of 
AGNs with z < 0.2 excluded. For the z > 0.2 sample, 
we could find good fits (with all of the KS probabilities 
in Lx7 z, and 2D exceeding 0.2) in any of the PLE and 
LDDE models. The acceptance of the PDE model was 
marginal (P2DKS ^ 0.1). The z < 0.2 regime is mainly 
contributed by AGNs in the RASS-based RBS and SA- 
N surveys, whose flux-area space have not been explored 
previously. Since these samples are completely identified 
(see Sect. ||) and we have included all emission-line AGNs, 
the relatively low value in this regime is not because of 
the incompleteness or sampling effects. The only source of 
possible systematic errors which could affect the analysis 
would be in the flux measurements, because of the differ- 
ences in details of the source detection methods among 
different samples. Some systematic shift of flux measure- 
ments might have occured between measurements in, e.g., 
the pointed and RASS data (for which there is no evi- 



dence). Thus we have made a sensitivity check by shift- 
ing the fluxes of all RBS and SA-N AGNs by -^20% and 
—20%. The flux-area relation (Fig. m has been modified 
accordingly. In either case in either value of fim, the ba- 
sic results did not change and especially the PLE model 
has been rejected with a large significance (with P2DKS 
ranging 10"^ — 10^^). 

4. Contribution to the Soft X-ray Background 

In this section, we discuss the contribution of AGNs to 
the soft X-ray background using the various models of the 
SXLF. As the absolute intensity level of the extragalac- 
tic 0.5-2 keV CXRB intensity, we use the results of an 
ASCA-ROSAT simultaneous analysis on the ASCA LSS 
field (Miyaji et al. in preparation), which covers a much 
larger field than Miyaji et al. ( 1998| ) and thus is subject 
to less uncertainties due to source fluctuations. There still 
are uncertainties in separation of the Galactic hard ther- 
mal and extragalactic components. Especially, it is still 
not clear whether the extragalactic component has also a 
soft excess at _E ^ 1 [keV] over the extrapolation from 
higher energies or whether the observed excess is domi- 
nated by the Galactic hard thermal component. Some au- 
thors prefer a model where the extragalactic component 
also contributes to the E '^1 [keV] excess because fit with 
a single power-law plus a thermal plasma would require 
an unusally low metal abundance of the thermal compo- 
nent for a Galactic plasma (Gendreau et al. |1995 ) and/or 
because many AGNs show soft excesses (e.g. Parmar et 
al. 1999). On the other hand, a self-consistent population 
synthesis model, including the AGN soft-excess below 1.3 
keV, still predicts that the low-energy excess is not promi- 
nent in the 0.5-2 keV range (Miyaji et al. 1999bD , mainly 
because the break energy shifts to the observed photon 
enrgy oi E ^ 0.4 [keV] for AGNs at z ~ 2, where the 
largest contribution to the CXRB is expected. The 0.25 
keV extragalactic component measured using a shadow- 



ing of a few nearby galaxies (Warwick & Roberts 1998) 



is consistent with both the single power-law extrapolation 
case and a slight soft excess (F ^ 2 for -E ^ 1 [keV]). 

In our comparison, we use (7.4 — 9.0) x 10^^^ 
[ergs^^ cm^^ deg~^] as a probable range of the extra- 
galactic 0.5-2 keV intensity, where the smaller value cor- 
responds to the single power-law form of the extragalactic 
component and the larger value corresponds to the case 
where the extragalactic component steepens to a photon 
index ofr = 2.3ati?^l [keV]. This range can be 
compared with the integrated intensity expected from the 
models. 

In Fig. O, we plot the cumulative soft X-ray (0.5-2 
[keV]) intensities of the model AGN populations as func- 
tions of redshift, /o.5-2kov(< z). As a reference, we have 
also plotted the cumulative contribution of the resolved 
AGNs in the sample, estimated by X^z <z '5'xiM('S'xi); 
where S'x, is the flux of the object i and A{Sxi) is the 
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Fig. 10. The cumulative 0.5-2 keV intensities /(< z) 
are plotted as a function of redshift for the PLE, PDE, 
LDDEl, and LDDE2 models for two different cosmologies 
as labeled. See caption for Fig. ^ for line styles corrcpond- 
ing to these four models. These curves include expected 
contribution from sources fainter than the survey limit 
using the model extrapolations. As a reference, the cumu- 
lative intensity /(< z) of the AGNs in the sample (see 
text) is also plotted (thin solid line with 90% errors) on 
each panel. This curve represents the contribution of ac- 
tually resolved and identified AGNs. Also the range of the 
0.5-2 keV extragalactic background intensity (see text) is 
shown by two horizontal thin dotted lines. 

available survey area at this flux (Fig. |]). The portion of 
the model curves above this line represents extrapolations 
to fainter fluxes than the limit of the deepest survey. 

It is apparent from Fig. |o| that the PDE model pro- 
duces almost 100% of the upper-estimate of the CXRB 
intensity, giving no room for, e.g. 10% contribution from 
clusters of galaxies (M99b), in the {n^,^,^) = (1.0,0.0) 
universe. In the low density universe with (f7ni,0A) = 
(0.3,0.0), the PDE model certainly overproduces the 
CXRB intensity. The PLE model produces about ~ 50% 
of the lower estimate of the CXRB in both cosmologies. 
The LDDEl model, which best describes the data in the 
observed regime, explains about 70% of the lower estimate 
of the CXRB intensity. The estimates are highly depen- 
dent on how one extrapolates the SXLF to fluxes fainter 
than the survey limit. In view of this, we explore an al- 
ternative LDDE model, which has been adjusted to make 
« 90% of the lower estimate of the extragalactic CXRB 
intensity, allowing ~ 10% contribution from clusters of 
galaxies. This version of the LDDE model (designated as 
LDDE2) has a fixed minimum evolution index Pmin in the 
LDDE formula. Then the first case of Eq. (|l^) is replaced 
by: 



e(z, Lx) = (1 + z)max(p„i„,pl-a(Log L,-Log L,)) ^ 



Table 4. Best-Fit LDDE2 Parameters 



{ilrn,^A) Parameters/KS probabilities 



LDDE2 
(1.0,0.0) 



LDDE2 
(0.3,0.0) 



LDDE2 
(0.3,0.7) 



;L. =0.85j:;g 



A = (0.88 ± .06) X 10" 

71 = 0.71 ± .14; 72 = 2.25 ± .10; pi = 4.82l:;^^ 

Zc = 1.64 ± .16;p2 = 0. (fixed) ;pmin = 4.0(fixed) 

Q = l.O(fixed); Log La = 44.1 (fixed) 

Pks = 0.7, 0.6, 0.3 (for L,z,2D) 



A 



.21 



(1.59 ± .10) X 10""; L, = 0.581 ^4 
71 = 0.55 ± .16; 72 = 2.30 ± .08; pi = 5.8 ± .3 
2c = 1.57±.12;p2 = 0.(fixed);pniin = 3.7(fixed) 
Q = 2.5 (fixed); Log La = 44.6 (fixed) 
Pks = 0-99, 0.5, 0.4 (for L,z,2D) 

A = (1.48 ± .09) X 10"^ L, = O.GOtil 

71 = 0.57 ± .17; 72 = 2.21 ± .08; pi = 5.3 ± .3 

2c = 1.59 ± .12;p2 = 0. (fixed) ;pmin = 3.3(fixed) 

Q = 2.5 (fixed); Log La = 44.5 (fixed) 

Pks = 0.7, 0.8, 0.4 (for L,2,2D) 



''Units- A: [/150 Mpc"^], L«: [lO''" /igo^erg s"^]. Parameter er- 
rors correspond to the 90% confindence level, search (see Sect. 

Oh. 



(2 < Zc;L^ < La). 



(11) 



We do not intend to represent a particular physical 
picture behind this formula. We rather intend to search 
for a formally simple expression which makes 90% of the 
CXRB and is still consistent with our sample in the regime 
it covers. We have searched for models accepted by the KS 
tests by adjusting parameters Pmin, a and La by hand and 
fitting by the maximum-likelihood method with respect to 
other variable parameters, requiring that the models give 
an integrated intensity of 6.7 x 10~^^ [erg s~^ cm~^ deg~^] . 
The parameter values of such LDDE2 models are listed in 
Table 4. 

By considering LDDE2, we have shown that there still 
is a resonable extrapolation of the AGN SXLF which 
makes up most of the soft CXRB. Of course this is not 
a unique solution. One may consider LDDEl and LDDE2 
as two possible extreme cases of how the SXLF can be 
extrapolated. Further implications are discussed in Sect. 



5. Evolution of luminous QSOs 

In this section, we consider QSOs with high soft X-ray 
luminosities (Log Lx > 44.5), where the behavior of the 
SXLF can be traced up to high redshifts. Also in the high- 
luminosity regime, at least in the local universe, we ob- 
serve very few absorbed AGNs, which could cause prob- 
lems with the K-correction, in the local universe (e.g. 
Miyaji et al. 1999b). If this tendency extends to the high 
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Fig. 11. The comoving number density of luminous 
(Log Lx > 44.5) QSOs in our ROSAT AGN sample are 
plotted as a function of redshift for two cosmologies as la- 
beled. The horizontal error bars indicate redshift bins and 
vertical error bars la errors. The top symbol of a down- 
ward arrow corresponds to the 90% (2.3 obj) upper limit. 
The points for (Om, ^a) = (0.3, 0.0) have been shifted hor- 
izontally by -1-0.1 in z for display purposes. The numbers 
of the X-ray luminnous QSOs for the four highest redshift 
bins are 24[32] (1.6 < z < 2.4), 8[12] (2.4 < z < 3.2), 
5[7] (3.2 < z < 4.6), and 0[0] (4.6 < z < 6.0) for 
(fim.fiA) = (1.0,0.0) [=(0.3,0.0)]. For comparison, num- 
ber density of optically-selected (Afe < —26) (dashed hue 
and filled tria ngles, from SSG95) and radio-selected (stars. 
Shaver et al. 1999|) QSOs, normalized to the soft X-ray 
selected QSO number density at z ^ 2.5 are overplotted. 
For the SSG95 data, this normalization corresponds to a 



multiplication by a factor of 7. Shaver et al. 1999 gave 
no absolute density.The optical and radio points are for 
(fi„„17A) = (1.0,0.0). 



redshift universe, our i?05'j4T' sample is a good representa- 
tion of luminous QSOs and the assumed single power-law 
of F = 2 would be a reasonable one. Thus we here investi- 
gate the evolution of the number density of the luminous 
QSOs using our sample. Fig. g shows that the SXLF at 
high luminosities can be approximated by a single power- 
law (7 « 2.3) for all redshift. Assuming this power-law 
with a fixed slope, we have calculated the number density 
of AGNs above this luminosity using the fitted normal- 
ization as described above in different redshift bins. The 
results are plotted in Fig. [Ill for two sets of cosmological 
parameters. Similar curves for optically and radio-selected 
QSOs are discussed below. 

In both cases, the number density increases up to 
z ~ 1.6 and flattens beyond this redshift. In both cos- 
mologies, the number density for z ^ 1.7 is consistent 
with no evolution. The Maximum-Likelihood fits in the 
z ^ 1.7, Log Lx > 44.5 region gave density evolution in- 



dices (ex [1 4- z]P) oi p = 0.5 ± 2.5 and p = 0.8 ± 2.1 
for (0„,0a) =(1.0,0.0) and (0.3,0.0) respectively. Subtle 
differences of the density curves seen in Fig. O between 
the two cosmologies come from two effects. Because differ- 
ent cosmologies give different luminosity distances, some 
objects which do not fall in the Log Lx > 44.5 region 
for (Oni,r2A) = (1.0,0.0) come into the sample in lower 
density cosmologies. Also the comoving volume per solid 
angle in a certain redshift range becomes larger in lower 
density cosmologies, thus the number density lowers ac- 
cordingly. These two effects work in the opposite sense 
and tend to compensate with each other, but the former 
effect is somewhat stronger. 

It is interesting to compare this curve with similar ones 
from surveys in other wavelengths. In Fig. |l^, we ov erplot 
number densities of optically- (Schmidt et al. 1995 , here- 
after SSG95) and radio-selected (Shaver et al. |l999D QSOs 
for (rim, ^a) — (1.0,0.0). The densities of these QSOs have 
been normalized to match the i?05'^T-selected QSOs at 
z '^ 2.5. This corresponds to a multiplicative factor of 7 



for the SSG95 sample. Shaver et al. (1999) gave no ab- 
solute number density. In order to assess the statistical 
significance of the apparent difference of the behavior at 
z > 2.7 between the ROSAT selected sample, we have 
used a Maximum-Likelihood fit to 17 QSOs in the sample 
in 2 > 2.2 and LogLx > 44.5. 



ex L^- ■ e{z) 



(12) 



(13) 



d$(£x,^) 
d Log Lx 

with 

( ■._ ( C (1.7<z<2.7) 

""^^^ ~\C exph/3(z - 2.7)] (2 > 2.7) 

where C is a constant. In above expression, (3 = corre- 
sponds to no evolution even for z > 2.7 and /3 = 1 is a good 
description of the rapid decrease of optically-selected QSO 
number density by SSG95. Fig. ^ shows that the radio- 
selected QSOs follow the SSG95 curve very well, but they 
do not have sufficient statistics to directly compare with 
the X-ray results. We have made a Maximum-Likelihood 
fit with only one free parameter: /3. Fixing 7 at 2.3, we 
have obtained the best-fit value and 90% errors (corre- 
sponding to A£ = 2.7) of /3 = 0.1^:5. The result changed 
very little if we treat 7 as a free parameter. Setting /? = 1 
increased the £ value by 3.3 from the best-fit value. This 
change in C corresponds to a 93% confidence level. The 
probability that (3 exceeds the value of 1 is « 4%, consid- 
ering only one side of the probability distribution. 

We have also checked statistical significance of the 
difference using the density evolution-weighted {V^/V^) 
statistics, (Avni & Bahcall 1980), which is a variant of 
the (l^'/T^^ax) statistics (Schmidt |196^ ) for the cases 
where surveys in different depths are combined. The 14 
and Vc are primed to represent that they are density 
evolution- weighted (comoving) volumes. If we take e(z) in 
Eq. {^ with /3 = 1 as the weighting function, (VJ/F^') 
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will give a value of 0.5 if the sample's redshift distri- 
bution follows the density evolution law of SSG95. An 
advantage of this method over the likelihood fitting is 
that one can check the consistency to an evolution law 
in a model-independent way, i.e., without assuming the 
shape of the luminosity function. Applying this statistics 
to 17 AGNs in z > 2.2, LogLx > 44.5, we have obtained 
(V^/V^) = 0.65 ± 0.07, where the la error has been es- 
timated by (12A^)~2 (A^:number of objects). The same 
sample has given unweighted (T4/T4) = 0.56 ± 0.07, con- 
sistent with a constant number density. If we use a harder 
photon index of F = 1.7 for the K-correction, 14 objects 
remain in this regime giving {V^/V^) = 0.60 ± 0.08. The 
inconsistency between the SSG95 optical results and our 
survey is thus marginal if high-redshift QSOs have a sys- 
tematically harder spectrum. 

6. Discussion 

In our analysis, we have found a good description of 
the behavior of the SXLF from a combination of various 
ROSAT surveys. As explained in Sect. 3.1, our expression 
is for the total AGN population, including type 1 and type 
2's and for the observed 0.5-2 keV band, because of the 
uncertainties of contents and evolution of AGNs in various 
spectral classes. These have to be assumed to find the best- 
bet K-corrected AGN evolution in the source rest frame. 
A detailed discussion of this aspect is beyond the scope 
of this paper. An approach for the problem is to make 
a population synthesis modeling, e.g., composed of unab- 
sorbed and absorbed AGNs similar to those of Madau et 



al. (|99J) and Comastri et al. ( |l995[ ) (see also Gilli |l999| for 
a recent work) . If one is constructing a model in a similar 
approach using our SXLF as a major constraint, what the 
model constructor should do is to calculate the expected 
SXLFs in the observed 0.5-2 keV band for all emission- 
line AGN populations (spectral classes) considered in the 
model (e.g. corresponding to different absorbing column 
densities) and then to compare the total model SXLF with 
our LDDE1/LDDE2 expressions. One version of our own 
models constructed using this approach has been shown 
in M99b. We do not recommend the use of the expressions 
in Appendix A. as the SXLF of unabsorbed AGNs for the 
reasons described there and Sect. 



3.1 



We have found two versions of LDDE expressions con- 
sistent with our sample in the luminosity and redshift 
regime covered: one which produces ~ 70% of the 0.5-2 
keV extragalactic CXRB (the lower etimate, see Sect. ^, 
and the other one which produces ~ 90%, as two rela- 
tively extreme cases on extrapolation. The real behavior 
is probably somewhere between these two. Note that we 
have only calculated the contribution to the CXRB for 
Log Lx > 41.7, where fits were made. Below this luminos- 
ity, we observe an excess (Fig. m, which conne cts we ll with 
the SXLF of near by Ga laxies (Schmidt et al. 1996| ; 



component has a local volume emissivity comparable or 
more to our sample AGNs in the 0.5-2 keV range and can 
contribute significantly to the soft CXRB. Because of the 
low luminosity, we can only detect this population in the 
very nearby universe in a large-area surveys like RASS. 
Deep small-area surveys would not give enough volume to 
detect them, since even the deepest part of RDS-LH can 
detect a Log Lx ~ 41 galaxy only up to z ~ 0.1. 

The X-ray emission of this low luminosity population 
is probably contributed by both star-formation and by 
low-activity AGNs (incl uding LINERS). Although Geor- 
gantopoulos et al. (1999)'s analysis suggests that a major 
contribution is from Seyfert galaxies and LINERS even 
at these low luminosities, star-formation activity can also 
contribute significantly to the X-ra y emis sion of these low- 
activity AGNs (see Lehmann et al. 1999a). As one extreme 
scenario, we assume that the X-ray emission from these 
low-luminosity sources is mostly from star-formation ac- 
tivity and their volume emissivity is assumed to evolve 
like the global star-formation rate (SFR; e.g. Madau et al. 



1996; Connolly et al. 1997), the integrated intensity would 



be roughly 30-40% of the lower estimate of the CXRB 
intensity. Even if the evolution of these low-luminosity 
sources were PLE, we would not detect any of them at 
intermediate to high redshifts even in the deepest ROSAT 
Survey on the Lockman Hole. Therefore this picture is 
still consistent with the result that the RDS-LH did not 
find any star burst galaxies. If the above scenario is the 
case, the bahavior of the AGN component would need 
to be close to LDDEl to allow room for a contribution 
from star-forming galaxies. In that case, the softer emis- 
sion from star-formation activity could contribute to the 
E ^ 1 [keV] excess of the CXRB spectrum and the total 
extragalactic 0.5-2 [keV] intensity could be closer to the 
upper estimate. If on the other hand, the large apparent 
local volume emissivity for the low-luminosity component 
is produced by the local overdensity and not representa- 
tive of the average present-epoch universe (e.g. Schmidt et 
al. 1996 is from a sample within 7.5 [Mpc]) and/or the X- 
ray evolution is slower than the global SFR (e.g. delayed 
formation of LMXB, White & Ghosh |1998| ), an LDDE2- 
like behavior for the Log L^ <; 41.7 AGN component may 
also be possible. A more detailed invetigation of the above 
scenarios and the exploration of other possibilities will be 
a topic of a future work. 

One of the most interesting results is the evolution of 
luminous QSOs discussed in Sect.^. A comparison of the 
evolution and the global star-formation rate is discussed 



Geor- 



gantopoulos et al. 1999| ), in the very local universe. This 



in Franceschini et al. (1999), where it is proposed that the 
evolution of the volume emissivity of the luminous QSOs 
evolves like the star-formation rate (SFR) of early-type 
galaxies, while that of the total AGN population (from 
the LDDEl and LDDE2 models) may evolve hke the SFR 
of all galaxies. Another interesting feature is that we find 
no evidence for a rapid decline of the QSO number den- 
sity at high redshift. The SSG95-like decrease at z > 2.7 is 
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marginally rejected. The difference may be caused by dif- 
ferent selection criteria. SSG95 have selected QSOs by the 
Lya luminosity and their QSOs are representative of more 
luminous QSOs (Mb < -26). Recently Wolf et al. ([1999D 



reported a similar tendency in their sample of QSOs from 
one of their CADIS fields, which typically have lower lumi- 
nosities than the SSG95 sample. Our X-ray selected AGNs 
with Log Lx > 44.5 have a seven times higher space den- 
sity than SSG95 at 2 ~ 2.5 and thus are sampling lower 
luminosity QSOs than SSG95. Thus if the behavior of our 
iZOS'^r-selected QSOs and those of Wolf et al. ( |l999D 
is really flat, this can be indicative of different formation 
epochs for lower and higher mass black holes. Adding more 
deep i? 05^1" surveys would enable us to trace the evolu- 
tion in this regime with a better statistical significance. 
The upcoming Chandra and XMM Surveys would extend 
the analysis to lower-luminsoity objects at the highest red- 
shifts as well as enabling us to give spectral information 
to separate the K-effect and the actual evolution of the 
number density. 

7. Conclusion 

We summarize the main conclusions of our analysis of the 
~ 690 AGNs from the i? 05^ T" surveys in a wide range of 
depths: 

1. Like previous works, we find a strong evolution of the 
SXLF up to z ~ 1.5 and a levelhng-off beyond this 
redshift. 

2. We have tried to find a simple analytical description 
of the overall SXLF. Our combined sample rejects the 
classical PLE model with high significance. The PDE 
model has been marginally rejected statistically and 
also overproduces the soft CXRB. 

3. We have found that an LDDE form (LDDEl), where 
the evolution rate is lower at low luminosities, gives an 
excellent fit to the overall SXLF. The extrapolation of 
the LDDEl form produces « 60 — 70% of the estimated 
extragalactic soft CXRB. 

4. Another form of LDDE (LDDE2), which equally well 
describes the overall SXLF from our sample, produces 
~ 90% of the extragalactic soft CXRB. These two 
LDDE models may be considered as two possible ex- 
treme cases when one considers the origin of the soft 
CXRB. 



greatly owes the effort of the ROSAT team and the optical fol- 
lowup teams in producing data and the catalogs used in the 
analysis. In particular, we thank K. Mason, A. Schwope, G. 
Zamorani, L Appenzeller, and I. McHardy for provoding us 
with and allowing us to use their data prior to publication of 
the catalogs. TM is supported by a fellowship from the Max- 
Planck-Society during his appointment at MPE. GH acknowl- 
edges DLR grants FKZ 50 OR 9403 5 and FKZ 50 OR 9908 
0. 



References 

Akiyama M., Ohta K., Yamada T. et al. 1999 in Proceeding 
of the first XMM workshop Science with XMM, in press 
( |astro-ph/9811012| ) 

Shanks T., Boyle B.J. et al. 1996, MNRAS 282, 



Almaini O. 

295 
Appenzeller I., Thiering I., Zickgraf F.-J. et al. 1998, ApJS 117, 

319 
Avni Y., Bahcall J.N. 1980, ApJ 235, 694 
Boldt E. 1987, Physics Reports 146, 215 
Bower R.G., Hasinger G., Castander F.J. et al. 1996, MNRAS 

281, 59 
Boyle B.J., Griffiths R.E., Shanks T., Stewart G.C., Georgan- 

topoulos I. 1993, MNRAS 260, 49 
Boyle B.J., Shanks T., Georgantopoulos I., Stewart G.C., Grif- 
fiths R.E. 1994, MNRAS 271, 639 
Comastri A., Setti G., Zamorani G., Hasinger G. 1995, A&A 

296, 1 
Connolly A., Szalay, A.S., Dickinson M., Subbarao M.U., Brun- 

ner, R.J. 1997, ApJ 486, Lll 
Delia Ceca R, Maccacaro T., Gioia I., Wolter A., Stocke J.T. 

1992 ApJ 389, 491 
Fasano G., Franceschini A. 1987, MNRAS 225, 155 
Fiore F., La Franca F., Giommi P. et al. 1999, MNRAS 306, 

L55 
Fischer J.-U., Hasinger G., Schwope A.D. et al. 1998, Astron. 

Nachr. 319, 347 
Franceschini A., Hasinger G., Miyaji T., Malquo li D. 1999, 



MNRAS (Letters), in press ( |astro-ph/9909290[ ) 
Gallagher S.C., Brandt W.N., Sam bruna R.M., Math ur S., Ya- 

masaki N. 1999, ApJ in press ( |:istro-ph/9902045| ) 
Gehrels N. 1986, ApJ, 303, 336 
Gendreau K.C., Mushotzky R.F., Fabian A.C. et al., 1995, 

PASJ 47, L5 
Georgantopoulos L, Basilakos S., Plionis M. 1999, MNRAS 



5. Th* evolution ot the number density ot luminous 
QSOs m our sample has been compared with that ot 
optically- and radio-selected QSOs. Our data are con- 
sistent with constant number density at z > 2.7, while 
optically- and radio-selected QSOs show a rapid de- 
cline. The statistical significance of this difference is 
just above 2a. Including more deep ROSAT surveys 
would trace the behavior with a better significance. 

Ackno wledgements. This works is based on a combination of 
extensl fe ROSAT surveys from a number of groups. Our work 



(Pink Pages) in press, ( |astro-ph/9902037| ) 
George I.M., Turner T.J., Netzer H. et al. 1998, ApJS 1 14, 73 
Gilh R., Risalti G. Salvati M. 1999 A&A in press, (|astro- 



ph/9904422) 



Griffiths R.E., Delia Cecca R., Georgantopoulos I. et al. 1996, 

MNRAS 281, 71 
Hasinger G., 1996 A&AS, 120, C607 
Hasinger G. 1998, Astron. Nachr. 319, 37 
Hasinger G., Burg R., Giacconi R., et al., 1993, A&A 275, 1 
Hasinger G., Burg R., Giacconi R., et al. 1998, A&A 329, 482 
Hasinger G., Lehmann L, Giacconi R. et al. 1999 in Highlights 

in X-ray Astronomy in Honor of Joachim Triimper's 65th 



Birthday, M PE Report (Garching:MPE) in press (astro- 



ph/9901103) 



Miyaji, Hasinger, Schmidt: 7?05'j4r AGN Luminosity Function 



17 



James F, 1994, MINUIT Reference Manual, CERN Program 

Library Long Writeup D506 (Geneva:CERN) 
Jones L.R., McHardy I.M., Merrifield M.R. et aL 1996, MN- 

RAS, 285, 547 
Lehmann L, Hasinger G., Schwope A.D., Boiler Th. 1999a, 

in Highlights in X-ray Astronomy in Honor of Joachim 

Trximper' s 65th Birthday, MPE Report (Garching:MPE), 

in press ( a,stro-ph/9810214 ) 
Lehmann L, Hasinger G., Schmidt M. et al. 1999b, A&A, sub- 
mitted 
Maccacaro T., Delia Cecca R., Gioia I.M. et al. 1991, ApJ 374, 

117 
Madau P., Ghisellini G., Fabian A.C. 1994, MNRAS 270, L17 
Madau P., Ferguson H.C., Dockinson M.E. et al. 1996, MNRAS 

283, 1388 
Marshall H.L., Avni Y., Tananbaum H., Zamorani G. 1983, 

ApJ 269, 35 
Mason K., Carrera F.J., Hasinger G. et al. 1999, MNRAS, 

submitted 
Mathur S., Elvis M., Singh K.P. 1995, ApJ 455, L9 
McHardy L, Jones L.R., Merrifield M.R. et al. 1998 MNRAS 

295, 641 
Miyaji T., Ishisaki Y., Ogasaka Y., Ueda Y. et al. 1998, A&A 

334, L13 
Miyaji T., Hasinger G., Schmidt M. 1999a, in Highlights in 

X-ray Astronomy in Honor of Joachim Trximper' s 65th 



Bir thday, M PE Report (Garching:MPE) in press, (astro 
ph/|9809398|, M99a) 



Miyaji T., Hasinger G., Schmidt M. 1999b, Adv. Sp. Res., in 

press (M99b) 
Newsam A.M., McHardy LM, Jones L.R., Mason K.O. 1998, 

MNRAS 292, 378 



astro 



Page M.J., Ca rrera F.J. 1999, MNRAS, in press ( 

ph/ p90943^ 
Page M.J., Carrera F.J., Hasinger G. et al. 1996, MNRAS 281, 

576 
Page M.J., Mason K.O., McHardy I.M., Jones L.R., Carrera 

F. J. 1997, MNRAS 291, 324 
Parmar A.N., Guainazzi, M., Oosterbroek T. et al. 1999 A&A, 

345, 611 
Piccinotti G., Mushotzky R.F, Boldt E.A. et al. 1982, ApJ, 

269, 423 
Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P. 

1992 Numerical Recipes in Fortran (Cambrideg: Cambridge 

Univ. Press), 640 
Romero-Colmenero E., Branduardi-Raymont G., Carrera F.J. 

et al. 1996, MNRAS 282, 94 
Schartel N., Schmidt M., Fink H. H. et al. 1997, A&A 320, 696 
Schmidt M. 1968, ApJ 151, 393 
Schmidt M, Green R.F. 1983, ApJ 269, 352 
Schmidt M., Schneider D.P., Gunn J.E. 1995, AJ 110, 68 

(SSG95) 
Schmidt K.-H., Boiler Th., Voges W. 1996, in 

Rontgenstrahlung from the Universe, MPE Report 

263 (Garching:MPE), 395 
Schmidt M., Hasinger G., Gunn J. et al. 1998, A&A 329, 495 
Schmidt M., Giacconi R., Hasinger G. et al. 1999, in X-ray 

Astronomy in Honor of Joachim Triimper's 65th Birthday, 

MPE Report, (Garching:MPE), in press 
Shaver P.A., Hook LM., Jackson C.A., WaU J.V., Kellermann 

K.L 1999, in Highly Redshifted Radio Lines eds. Carilh C. 

et al. (PASP: San Francisco), 163 



TuUy R.B., Shaya E.J. 1984, ApJ 281, 31 

Voges W. 1994, in Basic Space Science, Haubold H.J. & Onuora 

L.I. (eds.) 202 
Wisotzki L. 1998, Astr. Nachr. 319, 257 
White N. E., Ghosh P 1998, ApJ 504, L31 
Wolf C., Meisenheimer, K. Roser H.-J. et al. 1999, A&A 343, 

399 
Warwick R.S., Roberts T.P., 1998, Astron. Nachr. 319, 59 
Zamorani G., Mignoh M, Hasinger G. et al. 1999, A&A 346, 

731 
Zickgraf F.-J., Thiering I., Krautter J. et al. 1997, A&AS 123 

103 

Appendix A: SXLF of the 'type 1' AGN sample 

In the main part of this paper, we have concentrated on the 
SXLF expression for the mixture of type 1 and type 2 AGNs 
for the reasons explained in Sect. 3.1. However, since previous 
works in literature mainly give expression for only type 1 AGNs 
(with broad permitted lines) , it is of significant historical inter- 
est to investigate the SXLF properties for only type 1 AGNs. 
Because our samples come from several different sources and 
every subsample has its own criteria for classifying AGNs into 
subclasses, our expressions given here should not be used for 
any quantitative work (e.g. using it as a starting point of a 
population synthesis modeling under an assumption that they 
represent the unabsorbed AGNs) without assessment of possi- 



ble biases described in Sect, 3.1 



We have defined the 'type 1' AGN sample as follows. We 
have included AGNs explicitly classified in the original cata- 
logs as Seyfert 1-1. 5's, BLRGs, and QSOs, while excluded those 
classified as Seyfert 1.8-2, NELGs, and Narrow-line Seyfert I's 
(NLSl). The NLSl's have been excluded since they would not 
have been included in the 'broad-line' AGN samples in the 
previous works, especially those with low-quality optical spec- 
tra. A number of RBS objects classifed simply as 'AGNs' have 
been checked with the NED database and/or the original spec- 
tra for the subclassifications. For the Lockman Hole sample, we 
have included objects with ID classes (a)-(c) (see Schmidt et al. 
19981) and excluded (d)-(e). For t he M arano sample, those clas- 



sified as AGNs in Zamorani et al 1999 have been assumed to be 
type 1 AGNs unless otherwise stated, since type 2 AGNs have 
been explicitly noted. The AGNs which have not been sub- 
classified using the above procedure have been excluded from 
our 'type 1' sample. The fraction of AGNs included in this 
'type r sample are 98% (RBS), 90% (SA-N), 93%(RIXOS), 
88% (NEP-(-Marano-|-UKD), and 85% (LH) respectively. 

We have only considered the PLE and LDDEl models in 
this appendix. Table A.l shows best-fit para mete rs and KS 
probabilities (see main text for details). Table AT shows that 



PLE is still rejected with a large significance for both cosmolo- 
gies, while finding good fits with the LDDEl form. A plot sim- 
ilar to Fig.y is shown for the 'type 1' AGN sample in Fig. AT 
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Table A.l. Best-fit Parameters for the 'type 1' Sample 



Model 



PLE 

(1.0,0.0) 



PLE 

(0.3,0.0) 



LDDEl 

(1.0,0.0) 



LDDEl 

(0.3,0.0) 



Parameters/KS probabilities 



71 =0.43 ±.19; 72 = 2.30±.ll;pl = 3.0 ± .2 

2c = 1.45±.19;p2 = 0.3+g 

Pks = 5 10"*, 3 10-^,7 lO"'^ (for L,z,2D) 



A = (3.6 ± .2) X 10"^ L, = 0.34 ± .10 

71 =0.41 ±.19; 72 = 2.31±.ll;pl = 3.0 ± .2 

2c = 1.47 ± .28; p2 = 0.46 ± .7 

Pks = .02, .008, .002 (for L,2,2D) 

A = (1.40 ± .10) X 10""; L, = 0.60+2 

71 = 0.62 ± .20; 72 = 2.25 ± .09; pi = 5.4 ± .3 

2c = 1.55 ± .15; p2 = 0.0 (fixed) 

a = 2.5 ± .8; LogLa = 44.2 (fixed) 

Pks = 0.6,0.6,0.6 (for L,2,2D) ; 

A = (1.52 ± .10) X 10"''; L. = 0.55tJ^ 

71 =0.62 ±.23; 72 = 2.17 ±.08; pi = 5.3 ± .3 

2c = 1.62 ± .14; p2 = 0.0 (fixed) 

a = 3.0 ± .9; LogLa = 44.2 (fixed) 

Pks = 0.3,0.8,0.3 (for L,2,2D) 



See Captions for Tables H and H for units of the parameters 
and other notes. 
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Fig. A.l. Same as Fig. ^ except that this is for the 'type 
1' AGN sample and that only PLE and LDDEl models 
are plotted. 



